Heavy-residue isoscaling as a probe of the symmetry energy of hot fragments 
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The isoscaling properties of isotopically resolved projectile residues from peripheral collisions of 
*'^Kr (25 MeV/nucleon), '^''Ni (25 MeV/nucleon) and "''Xe (20 MeV/nucleon) beams on various 
target pairs are employed to probe the symmetry energy coefficient of the nuclear binding energy. 
The present study focuses on heavy projectile fragments produced in peripheral and semiperipheral 
collisions near the onset of multifragment emission (E*/A = 2-3 MeV). For these fragments, the 
measured average velocities are used to extract excitation energies. The excitation energies, in turn, 
are used to estimate the temperatures of the fragmenting quasiprojectiles in the framework the 
Fermi gas model. The isoscaling analysis of the fragment yields provided the isoscaling parameters 
a which, in combination with temperatures and isospin asymmetries provided the symmetry energy 
coefficient of the nuclear binding energy of the hot fragmenting quasiprojectiles. The extracted 
values of the symmetry energy coefficient at this excitation energy range (2-3 MeV/nucleon) are 
lower than the typical liquid-drop model value --^ 25 MeV corresponding to ground-state nuclei and 
show a monotonic decrease with increasing excitation energy. This result is of importance in the 
formation of hot nuclei in heavy-ion reactions and in hot stellar environments such as supernova. 

PACS numbers: 25.70.Mn,25.70.Lm,25.70.Pq 



I. INTRODUCTION 

The study of the nuclear symmetry energy is currently 
a topic of intense theoretical and experimental work. It is 
well established that the symmetry energy plays a central 
role in a variety of astrophysical phenomena including the 
structure and evolution of neutron stars and the dynam- 
ics of supernova explosions In 
addition, the symmetry energy determines the nuclear 
structure of neutron-rich or neutron deficient rare iso- 
topes HIISIIllil- 

The symmetry energy at normal nuclear density has 
been obtained from a number of many-body approaches 

in III in HI 113 

and from nuclear mass systematics 
[isL ligl V2Sx l2ll |. However, its values at densities below 
or above the normal nuclear density are not adequately 
constrained j22j 2^ 2M ■ Indeed, the experimental deter- 
mination of the symmetry energy and its density depen- 
dence is a challenging scientific endeavor. Information 
on the symmetry energy can be gleaned from the deter- 
mination of the neutron skins of neutron-rich nuclei 0, 
from elastic scattering on neutron-rich nuclei and 
from heavy-ion collisions. For the latter, a great deal of 
effort is currently devoted in order to identify observables 
sensitive to the nuclear symmetry energy and its density 
dependence (see, e.g., lUlillHlll'ZSj). 

One important observable in heavy ion collisions is the 
fragment isotopic composition investigated with the re- 
cently developed isoscaling approach 30] . The isoscaling 
approach attempts to isolate the effects of the nuclear 
symmetry energy in the fragment yields, thus allowing a 



direct study of the role of this term of the nuclear bind- 
ing energy in the formation of hot fragments. Isoscaling 
refers to a general exponential relation between the yields 
of a given fragment from two reactions that differ only in 
their isospin asymmetry (N/Z) [Hl^ls^- In particu- 
lar, if two reactions, 1 and 2, lead to primary fragments 
having approximately the same temperature but differ- 
ent isospin asymmetry, the ratio R2i(N,Z) of the yields of 
a given fragment (N,Z) from these primary fragments ex- 
hibits an exponential dependence on the neutron number 
N and the atomic number Z of the form: 



R2i{N, Z) = Cexp{aN + PZ) 



(1) 
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where a and P are the scaling parameters and C is a nor- 
malization constant. This scaling behavior has been ob- 
served in a very broad range of reactions including evap- 
oration J33, , fission (3^ , deep- inelastic reactions 
JMJsl Hg and multifragmentation [M HI IHIll El 
[ifllU. In the initial studies of isoscaling, it was shown 
that the isoscaling parameters are almost unaffected by 
the sequential decay of the primary fragments, due to 
possibly similar de-excitation paths of the two primary 
fragments, thus they could provide information on the 
early stage of fragment formation. In particular, within 
the statistical framework, the isoscaling parameter a is 
linearly related to the symmetry ener gy c oefficient of the 
fragment binding energy |2l El ElliS ■ This relation 
provides the key connection of the measured isoscaling 
parameter with the symmetry energy coefficient. 

In the present work, the isoscaling properties of isotopi- 
cally resolved projectile residues from peripheral colli- 
sions of^^Kr (25 MeV/nuclcon), (25 MeV/nucleon) 
and i^^Xe (20 MeV/nucleon) beams on a variety of tar- 
get pairs are employed to probe the symmetry energy 
coefficient of the nuclear binding energy. The collection 
and complete characterization of the residues in terms 
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of their atomic number Z, mass number A and veloc- 
ity has been performed with the use of two magnetic 
separators: the MARS recoil separator and the Super- 
conducting Solenoid Line at Texas A&M University. In 
this work, apart from isotopically resolved yields, the 
velocities of the fragments are obtained with high res- 
olution and are used to provide information on the ex- 
citation energy (and temperature) of the primary frag- 
ments. The isoscaling parameters a along with temper- 
atures and isospin asymmetries yielded the values of the 
nuclear symmetry energy in the excitation energy range 
2-3 MeV/nucleon. The paper is organized as follows. 
In Section II, a brief description of the experimental de- 
vices, the measurements and the data analysis is given. 
In Section III, the isotopic scaling of the fragment yields 
and the velocity distributions are presented. In Section 
IV, the systematics of the isoscaling parameter a with 
respect to isospin asymmetry is presented and used to 
get the symmetry energy. Finally, conclusions from the 
present study are summarized in Section V. 



II. EXPERIMENTAL METHODS AND DATA 
ANALYSIS 

The present studies were performed at the Cyclotron 
Institute of Texas A&M University using two different 
devices: the MARS recoil separator and the Supercon- 
ducting Solenoid Line. Below we give a concise descrip- 
tion of the measurements with each of these devices. 



A. Measurements with the MARS recoil separator 

The reactions of ^^Kr with ^^'^^Ni and I24,ii2sn were 
studied with the MARS recoil separator. The general 
isoscaling analysis of these data has been reported re- 
cently [33 • Nevertheless, for a complete presentation of 
the heavy-residue isoscaling approach to probe the sym- 
metry energy, we briefly summarize the method below. 
A 25 MeV/nucleon *^Kr beam from the K500 supercon- 
ducting cyclotron, with a current of ^^1 pnA, interacted 
with isotopically enriched targets of ^^Ni, ^^Ni and ^^''Sn, 
^^^Sn. The reaction products entered the MARS spec- 
trometer |43| having an angular acceptance of 9 msr and 
momentum acceptance of 4%. The primary beam struck 
the target at 4.0° relative to the optical axis of the spec- 
trometer. Fragments were accepted in the polar angular 
range 2.7°-5.4°. This angular range lies inside the graz- 
ing angle 9gr=6-5° of the Kr-|-Sn reactions and mostly 
outside the grazing angle ^^-=3.5° of the Kr-fNi reac- 
tions at 25 MeV/nucleon >44|. (The spectrometer angle 
setting was chosen to optimize the production of very 
neutron-rich fragments from the Kr-|-Sn systems whose 
detailed study has been reported in 0-) ^1 strip- 
per foil (1 mg/cm^) was used to reset to equilibrium the 
ionic charge states of the projectile fragments. MARS op- 
tics [4^ provides one intermediate dispersive image and 



one (final) achromatic image (focal plane). At the focal 
plane, the fragments were collected in a 5x5 cm two- 
element (AE, E) Si detector telescope. Time of flight was 
measured between two PPACs (Parallel Plate Avalanche 
Counters) positioned at the dispersive image and at the 
focal plane, respectively, and separated by a distance of 
13.2 m. The PPAC at the dispersive image was also X- 
Y position sensitive and used to record the position of 
the fragments. The horizontal position, along with NMR 
measurements of the field of the MARS first dipole, pro- 
vided the magnetic rigidity. Bp, of the particles. 

The determination of the atomic number Z was based 
on the energy loss of the particles in the AE detec- 
tor and their velocity. The ionic charge q of the parti- 
cles entering MARS was obtained from the total energy 
Etot=AE -I- E, the velocity and the magnetic rigidity. 
The measurements of Z and q had resolutions of 0.5 and 
0.4 units (FWHM), respectively. Since the ionic charge 
is an integer, we assigned integer values of q for each 
event by putting windows {Aq — 0.4) on each peak of 
the q spectrum. Using the magnetic rigidity and veloc- 
ity measurement, the mass-to-charge A/q ratio of each 
ion was obtained with a resolution of 0.3%. Combin- 
ing the q determination with the A/q measurement, the 
mass A was obtained as: A — qint x ^/g {(lint is the 
integer ionic charge) with a resolution (FWHM) of 0.6 
A units. Combination and normalization of the data 
at the various magnetic rigidity settings of the spec- 
trometer (in the range 1.3-2.0 Tm), summation over 
all ionic charge states and, finally, normalization for 
beam current and target thickness, provided fragment 
yield distributions with respect to Z, A and velocity. 
Further details of the analysis procedure can be found 
in j4^ (4^. The yield distributions, summed over ve- 
locities, were used to obtain the fragment yield ratios 
R2i{N, Z) = Y2{N, Z)/Yi{N, Z) employed in the present 
isoscaling studies. 



B. Measurements with the Superconducting 
Solenoid Line 

The heavy-residue work with the Superconducting 
Solenoid Line shares many similarities with the previ- 
ously described work with MARS and it is described be- 
low. The complete (two-stage) Superconducting Solenoid 
Line (Fig. 1) consists of the 7-Tesla superconducting 
solenoid ( "BigSol" ) of the University of Michigan (first 
stage) "is, 49j and a large-bore quadrupole triplet (sec- 
ond stage). The whole separator line is also referred to 
as the "BigSol" line. Details of the development of the 
line and plans/progress towards producing rare isotopes 
in deep- inelastic collisions are provided in 47, 50.] . 

In the present isoscaling studies with BigSol, first, the 
reactions of a ^^Ni (25 MeV/nucleon) beam with targets 
of 64,58Nj^ 124,ii2Sjj 208pb^ 232rpj^ ^^^.g gtudied. Sec- 
ond, the reactions of a ^'^^Xe (20 MeV/nucleon) beam 
with targets of ^^^^^Wi, I24,ii2sn ^nd ^^^Au, 232Th were 
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were performed. The data were normalized and appro- 
priately combined, following the procedure described in 
pq . Summation over all ionic charge states and, finally, 
normalization for beam current and target thickness, pro- 
vided the fragment yield distributions with respect to Z, 
A and velocity. As in the case of the MARS data, the 
fragment yield distributions were summed over velocities 
and used to obtain the yield ratios for isoscaling. 



FIG. 1: Schematic diagram of the Superconducting Solenoid 
Line (BigSol) used in the present heavy residue isoscaling 
studies (see text) 

measured. The typical beam current was '--^l pnA. The 
primary beam, after hitting the target, was collected on 
an on-axis blocker located ^30 cm after the target. The 
beam blocker along with a circular aperture at this lo- 
cation defined the angular acceptance of the line to be 
1.5°-3.0°. The fragments then passed through a stripper 
foil, traversed the solenoid and were focused at the in- 
termediate image (Fig. 1). At this location, a magnetic 
rigidity (or momentum-over-charge p/q) acceptance of 
~4% was defined with another circular aperture. Sub- 
sequently, the fragments were transported through a 7.5 
meter line and focused with the aid of the quadrupole 
triplet at the end of the line (final image) . Time-of-flight 
was provided between two X-Y position sensitive PPACs, 
one at the intermediate image and the other at the fi- 
nal image. A silicon detector array similar to the one 
used in the MARS measurements provided AE, E infor- 
mation, which, combined with time-of-flight (with res- 
olution ~0.5%) provided Z and A determination (with 
resolutions of 0.5 and 0.6 units, respectively, for Ni-like 
fragments) . 

It should be noted that, in contrast to the MARS mea- 
surements, here the mass determination was based solely 
on total energy and time-of-flight. The reason is that Big- 
Sol does not provide high- resolution p/q dispersion. The 
Solenoid Line is not a magnetic spectrometer in the clas- 
sical sense Its elements are only focusing elements. 
There is radial p/q dispersion at the intermediate image 
due to the solenoid (and, of course, at the final image, 
due to the quadrupole triplet) as a consequence of the 
variation of the location of the focus with p/q. This dis- 
persion is also a function of the initial angle The 
measurement of the radial distance of the fragments at 
the intermediate image (emerging, as stated previously, 
in the initial angular range of 1.5°-3.0°), combined with 
the value of the central magnetic field of the solenoid, 
provided a p/q determination with a resolution of ^^2%. 
Even though this resolution is not useful to improve the 
mass determination (as was done in the MARS data [4fil|). 
it was adequate to specify the charge state q of the ions. 

As in the case of the MARS data, a series of runs at 
overlapping magnetic rigidity settings of the line in the 
range 1.1-1.6 Tm for the ^''Ni (25MeV/nucleon) data 
and 1.0-1.5 Tm for the ^^^Xe (20 MeV/nucleon) data 



III. ISOSCALING AND EXCITATION ENERGY 
DATA FROM THE VARIOUS REACTIONS 

From the measured fragment yield distributions, we 
construct the ratio R2i{N,Z) = Y2{N, Z)/Yi{N, Z) of 
yields of a given projectile fragment (N,Z) from reactions 
2 and 1 using the convention that index 2 refers to the 
more neutron-rich system and index 1 to the less neutron- 
rich one. The results of the reactions with each of the 
various beams are given in the following subsections. 



A. '^''Kr (25MeV/nucleon) data (MARS) 

Fig. 2 shows the yield ratios R2i(N,Z) as a function of 
fragment neutron number N for selected isotopes for the 
Kr-|-Ni reactions (top panel) and the Kr-t-Sn reactions 
(bottom panel). The different isotopes are shown by al- 
ternating filled and open symbols for clarity. For each 
element, an exponential function of the form Gexp{a N) 
was fitted to the data and is shown in Fig. 2 for the 
selected elements. In the semi-log representation, the 
straight lines for each element are parallel up to Zr^34 
for Kr-|-Ni and up to Z~28 for Kr-|-Sn. For heavier frag- 
ments from the Kr-|-Sn system, the fits to the data show 
gradual decrease in the slopes with increasing Z, a be- 
havior that has been shown to manifest incomplete N/Z 
equilibration for the most peripheral events for this reac- 
tion whose projectile-like fragments were observed inside 
the grazing angle [35ll36| . 

In Fig. 3a, we present the slope parameters a of the 
exponential fits (obtained as described for Fig. 2) as 
a function of Z. For the Kr-|-Ni reactions (open sym- 
bols), the slope parameter a is constant in the whole 
range Z=12-36 at an average value a=0.254. For the 
Kr-|-Sn reactions, the parameter a is roughly constant 
with an average value of 0.415 in the range Z=12-28 and 
then it decreases for Z>28. In the case of the Kr-fNi 
systems under the present experimental conditions (ob- 
servation outside the grazing angle), isoscaling holds es- 
sentially in the whole range of observed fragments. To 
obtain information about the excitation energy of the 
primary fragments from the present reactions, we will 
employ the correlation of the measured velocity with the 
atomic number [3^. Fig. 3b presents the average veloc- 
ities of the fragments as a function of Z. Open symbols 
correspond to the ®^Kr-|-^'*Ni reaction and closed symbols 
to the ^^Kr-f ^^''Sn reaction. In this figure, we observe 
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FIG. 2: (MARS data) Yield ratios R2i{N,Z) = 
Y2(N, Z)/Yi{N, Z) of projectile residues from the reactions 
of *®Kr (25 MeV/nucleon) with ^^'^^Ni (a), and ^^^'"^Sn (b) 
with respect to N for the Z's indicated. The data are given 
by alternating filled and open circles, whereas the lines are 
exponential fits. 



that for fragments close to the projectile, the velocities 
are slightly below that of the projectile, corresponding to 
very peripheral, low-excitation energy events. A mono- 
tonic decrease of velocity with decreasing Z is observed, 
indicative of lower impact parameters, and thus, higher 
excitation energies. For the ^^Kr+^^'^Sn reaction (closed 
symbols in Fig. 3b), the descending velocity-Z correla- 
tion continues down to Z^28; for lower Z's, the velocity 
starts increasing with decreasing Z. A minimum veloc- 
ity for Z'~28 can be understood by assuming that these 
residues originate from primary fragments with a maxi- 
mum observed excitation energy. Fragments with Z near 
the projectile down to Z~28 originate from evaporative 
type of deexcitation which does not modify, on average, 
the emission direction of the residues. Thus, the residue 
velocities provide information on the excitation energy. 
Residues with lower Z arise from primary fragments un- 
dergoing cluster emission and/or multifragmentation and 
the velocity of the inclusively measured fragments is not 
monotonically related to the excitation energy. For the 
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FIG. 3: (MARS data) (a) Isoscaling parameter a as a func- 
tion of Z for projectile residues from the reactions of '^''Kr 
(25 MeV/nucleon) with ^'"^'•''^Ni (open circles) and ^^'''"^Sn 
(closed circles). The straight lines are constant value fits 
for each system, (b) Average velocity versus atomic num- 
ber Z correlations for projectile residues from the reactions of 
^''Kr (25 MeV/nucleon) with "^^Ni (open symbols) and ^^^Sn 
(closed symbols). The dashed line (marked "PR") gives the 
velocity of the projectile, whereas the arrows indicate the min- 
imum average residue velocities observed. 



^®Kr-|-^''Ni reaction, a similar behavior is observed. How- 
ever, the decreasing velocity-Z correlation is observed 
down to Z~32. We remind that fragments from this re- 
action were measured mostly outside the grazing angle, 
so that they correspond to more damped collisions, in 
such a way that the final residues receive a larger recoil 
during the deexcitation and appear within this angular 
range. For Z~30 32, we observe a minimum velocity and 
for lower Z's an increase of the velocity with decreasing 
Z analogous to the ^^Kr-|-^^^Sn reaction. We mention 
that the average velocities from the reactions with the 
neutron-poor targets are, within the experimental uncer- 
tainties, almost the same as the corresponding from the 
reactions with the neutron-rich targets for both pairs of 
systems and are not shown in Fig. 3b. For both reac- 
tions, the ascending part of the velocity vs Z correlation 
for the lower part of the Z range is primarily due to the 
combined effect of angle and magnetic rigidity selection 
imposed by the spectrometer. 

Employing the observed minimum velocities for the 
Kr-|-Ni and Kr-|-Sn reactions and, furthermore, assum- 
ing two-body kinematics, we can estimate the total ex- 
citation energy of the quasiprojectile - quasitarget sys- 
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tern employing standard mass tables |53|- In regards to 
the sharing of excitation energy, a reasonable assumption 
for per i p heral /semiperipheral collisions is equal division 
|53>IS3>IS3 at relatively low kinetic energy losses, Ejo^s- 
This assumption was employed in our previous analysis 
of these isoscaling data 35]. In the present work, we use 
a more appropriate prescription in agreement with the 
experimentally observed transition of the excitation en- 
ergy sharing from the equal division limit (at low E;oss) 
to the thermal limit that may be attained near the max- 
imum of the kinetic energy loss, i^ioss,max- Specifically, 
we estimate the fraction of the excitation energy of the 
quasi-projectile assuming a linear evolution, with respect 
to T^ioss/^ioss.max, from the equal division limit to the 
thermal equilibrium limit. Following the above proce- 
dure, we can estimate an average excitation energy per 
nucleon for the hot quasiprojectile fragments of E*/A = 
2.4 MeV for the Kr+Ni system and E*/A = 2.0 MeV for 
the Kr-|-Sn system (Table I). We note that, for these sys- 
tems, the present estimates are close to the value E*/A 
= 2.2 MeV obtained under the assumption of equal ex- 
citation energy division employed in our previous work 



B. ^-"Ni (25MeV/nucleon) data (BigSol) 

In Figs. 4 and 5, we show the data for the ^''Ni (25 
MeV/nucleon) reactions obtained with the BigSol line. 
The presentation of the data follows a line similar to 
that of the MARS data (Figs. 2 and 3). In Figs. 4a,b, 
a general isoscaling behavior is seen in both the Ni-|-Ni 
(grazing angle 6*3^=3.8°) and Ni-|-Sn (0^^=6.5°) reactions 
essentially in the whole range of elements measured in 
the magnetic rigidity range 1.1-1.6 Tm. A slight de- 
crease of the isoscaling parameter a is discernible for 
near-projectile elements (Fig. 5a) possibly due to in- 
complete N/Z equilibration. Constant value fits to the a 
parameter vs Z in the range Z=12-24 yielded a=0.250 
and 0.324 for Ni+Ni and Ni-hSn respectively (Table I). 

For the ^-^Ni + '^^'^Th, ^ospb pair of reactions 
{9gr=9.5°), we also observe isoscaling despite the small 
difference of ~2.5% in the N/Z of the two systems (com- 
pared to, e.g., the ~10% N/Z difference in the reaction 
pair of ^''Ni -I- ^^'*Sn,^^^Sn). This limiting case shows the 
sensitivity of the isoscaling signal to the N/Z of the re- 
acting systems. A constant value fit to the a parameter 
data of the Ni-|-Th,Pb pair in the region Z=14-18 (Fig. 
5a) gives, as expected, a small value a=0.060 (Table I). 

The average velocities of the three ^''Ni reactions are 
shown in Fig. 5b (only the neutron-rich systems are pre- 
sented). As we see in the figure, near-projectile residues 
with velocities close to that of the projectile were not 
collected in these measurements, mainly due to magnetic 
rigidity selection (and, in addition, due to angle selec- 
tion A6'=l. 5-3.0°, particularly for the Ni+Th,Pb sys- 
tems with 9gr=9.5°). In a manner similar to the MARS 
data, we employed the observed minimum velocities to 
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FIG. 4: (BigSol data) Yield ratios R2i{N, Z) = 
Y2(N , Z) /Y\{N , Z) of projectile residues from the reactions 
of '^''Ni (25 MeV/nucleon) with ^^'^^m (a), ^^"^'^^Sn (b), and 
232rpjj 208p^j 1^^^ ^j^jj respect to N for the Z's indicated. The 
data are given by alternating filled and open circles, whereas 
the lines are exponential fits. 



extract average excitation energies for the fragmenting 
Ni-like quasiprojectiles resulting in a common value of 
E7A = 2.9 MeV for the three pairs of the ^''Ni reactions 
(Table I). 



C. "''Xe (20 MeV/nucleon) data (BigSol) 

Finally, in Figs. 6 and 7, we present the data for the 
reactions of the ^'^^Xe (20 MeV/nucleon) beam on the 
three pairs of targets 64.58Ni, 1244123^ ^nd 232Th, ^^Uu 
obtained with the BigSol line. Again, the presentation 
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FIG. 5: (BigSol data) (a) Isoscaling parameter a as a func- 
tion of Z for projectile residues from the reactions of ^''Ni (25 
MeV/nucleon) with ^^-''^Ni (open circles), ^^^'"^Sn (closed 
circles), and 232Th,2°*Pb (open diamonds). The straight 
lines are constant value fits for each system, (b) Average 
velocity versus atomic number Z correlations for projectile 
residues from the reactions of ^'*Ni(25 MeV/nucleon) with 
'^"Ni (open symbols), ^^''Sn (closed symbols), and ^sa-j-jj (open 
diamonds) . The dashed line (marked "PR" ) gives the velocity 
of the projectile. 



of the data is similar to that of the preceding sections. 
Fig. 6 shows the isoscahng behavior of the yields and 
the fits to selected elements. Fig. 7a shows the variation 
of the isoscaling parameter a as a function of Z. The 
a values for Xc+Ni and Xc+Sn arc displaced vertically 
for viewing. Despite the relatively large fluctuation in 
the data points, the data were fitted with straight lines 
in the region Z=14 46 yielding the values: 0.129, 0.193, 
and 0.096 for the three reaction pairs, respectively (Table 
I)- 

Fig. 7b shows the average velocities of the three -'^^^Xe 
reactions (again, only the neutron-rich systems are pre- 
sented). The grazing angles for these systems are 3.9°, 
6.3° and 10.0° respectively. The angle and magnetic 
rigidity selection were such that near-projectile residues 
with velocities close to that of the projectile; were not 
collected, as in the case of the ^^Ni reactions. However, 
the kinematical conditions were such that lighter frag- 
ments (below Z=35) show an ascending velocity behav- 
ior with decreasing Z analogous to the one observed in 
the Kr-|-Ni,Sn reactions (Fig. 3b). Using the observed 
minimum velocities, average excitation energies for the 
fragmenting Xe-like quasiprojectiles were extracted re- 
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(a) 136Xe + 64,58Ni 
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(b) i36xe + i24,ii2sn 
Z = 14 18 22 26 30 34 38 42 46 
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(c) i36xe + 232Th, is^Au 

Z = 14 18 22 26 30 34 38 42 46 
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FIG. 6: (BigSol data) Yield ratios R2i{N,Z) = 
Y2{N, Z)/Yi{N, Z) of projectile residues from the reactions of 
i36Xe (20 MeV/nucleon) with "^^-^^Ni (a), i24.ii2sn (b), and 
232rj,jj 197^^ (c) with respect to N for the Z's indicated. The 
data arc given by alternating filled and open circles, whereas 
the lines are exponential fits. 



suiting again in a common value of E*/A = 2.5 MeV for 
the three pairs of the ^^^Xe reactions. 



IV. INTERPRETATION OF RESULTS AND 
DISCUSSION 

Before proceeding to the analysis and interpretation 
of the present results on heavy fragment isoscaling, we 
present a summary of the reactions and the relevant pa- 
rameters in Table I. The difference in the isotopic compo- 
sition, expressed as the proton fraction squared, (Z/A)"^, 



6Xe (20MeV/u) + Ni, Sn, Th-Au 



0.80 
0.70 



0.50 - ^tOT I * 



0.40 
0.30 
0.20 



(a) 
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FIG. 7: (BigSol data) (a) Isoscaling parameter a as a function 
of Z for projectile residues from the reactions of ^^®Xe (20 
MeV/nucleon) with '^*''^^Ni (open circles), ^'^^'^^'^Sn (closed 
circles), and ^^^Th,^''**Pb (closed diamonds). The values of 
the first two systems are displaced vertically by +0.4 and +0.2 
units, respectively, for viewing. The straight lines are con- 
stant value fits for each system, (b) Average velocity versus 
atomic number Z correlations for projectile residues from the 
reactions of ^^"^Xe (20 MeV/nucleon) with '^''Ni (open sym- 
bols), ^^''Sn (closed symbols), and ^'■'^Th (open diamonds). 
The dashed line (marked "PR") gives the velocity of the pro- 
jectile. 



for each pair of systems is given along with the excita- 
tion energies (obtained from residue velocities), as well as 
the temperatures (calculated with the Fermi gas model 
and the expanding mononucleus model , as discussed 
below). The isoscaling parameters a and, finally, the 
extracted values of the symmetry energy coefficient are 
summarized. 



TABLE I; Summary of reaction pairs studied in this work 
(see text), along with parameters relevant to the isoscaling 
analysis: A{Z/A)'^ difference in (Z/A)^ of the two reactions 
of each pair. E*/A (MeV): excitation energy per nucleon 
of the corresponding primary quasiprojectiles. K (MeV): in- 
verse level density parameter. Tf, Tm(MeV): temperature 
estimates obtained from the Fermi gas model and the ex- 
panding mononucleus model [s^ . a: isoscaling parameter, 
c = a/A{Z/A)^ : reduced isoscaling parameter. Csym,,F, 
Csym,m (MeV): Symmetry energy coefficient of the nuclear 
binding energy obtained using the two temperature determi- 
nations Tf and T™ respectively, as mentioned above. Errors 
(one standard deviation) of the measured quantities are given 
in parentheses. 



Reaction 


A(Z/A)2 


E*/A 


K 


Tf 




a c 




Csym,m 


Kr+Sn 


0.0209 


2.0 


11.5 


4.8 


3.9 


0.415 19.9 


23.9 


19.4 






(0.1) 








(0.010) (0.5) 


(0.6) 


(0.5) 


Kr+Ni 


0.0154 


2.4 


11.8 


5.3 


4.2 


0.254 16.5 


21.9 


17.2 






(0.1) 








(0.005) (0.3) 


(0.5) 


(0.4) 


Ni+Ni 


0.0193 


2.9 


11.7 


5.8 


4.3 


0.250 13.0 


18.9 


14.0 






(0.1) 








(0.007) (0.4) 


(0.6) 


(0.4) 


Ni+Sn 


0.0240 


2.9 


11.7 


5.8 


4.3 


0.324 13.5 


19.6 


14.5 






(0.1) 








(0.005) (0.3) 


(0.5) 


(0.4) 


Ni+Th-Pb 


0.0046 


2.9 


11.7 


5.8 


4.3 


0.060 13.0 


18.9 


14.0 






(0.1) 








(0.010) (2.2) 


(3.2) 


(2.2) 


Xe+Ni 


0.0106 


2.5 


13.5 


5.8 


4.5 


0.129 12.2 


17.7 


13.7 






(0.1) 








(0.006) (0.6) 


(0.9) 


(0.7) 


Xe+Sn 


0.0159 


2.5 


13.5 


5.8 


4.5 


0.193 12.1 


17.6 


13.6 






(0.1) 








(0.005) (0.3) 


(0.5) 


(0.4) 


Xe+Th-Au 


0.0064 


2.5 


13.5 


5.8 


4.5 


0.096 15.0 


21.8 


16.9 






(0.1) 








(0.005) (0.8) 


(1.2) 


(0.9) 



A. Determination of temperature 

In order to estimate the temperature of the fragment- 
ing quasiprojectiles using the measured average excita- 
tion energies, we first employ the simple Fermi gas re- 
lationship of the form E* — ^T"^, with T the tempre- 
tature and K the inverse level density parameter. It is 
well known that for the non-interacting Fermi gas model 
the value of the inverse level density parameter is K ~16 
(e.g. [131), whereas the experimental data are consistent 
with lower values of K dependent on both the excitation 
energy and the mass range |53 |. For the present sys- 



tems in the 2-3 MeV/nucleon excitation range, values in 
the vicinity of K=12-13 are consistent with the exper- 
imental systematics '58]. Given the expected mass and 
excitation energy dependence of K, we decided to use 
values of K following the model of Shlomo and Natowitz 
[5^ |H3| . This model determines the nuclear level density 
within the framework of the Fermi gas incorporating the 
effects of the finite size of the nucleus, the contributions 
of the continuum states and the temperature and den- 
sity dependence of the nucleon effective mass. For the 
hot quasiprojectiles (Kr-like, Ni-like and Xe-like) of the 
present work, the values of K reported in Table I were 
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obtained by interpolating the results presented in Fig. 
1 of |52li • The corresponding values of the temperature 
are also listed under the column and discussed later 
in relation to Figs. 9b and 9c. These temperature val- 
ues are in reasonable overall agreement with measured 
temperatures of similar systems in this excitation e nerg y 
range, as systematically analyzed and presented in |58l |. 

In addition to the above procedure to obtain the tem- 
perature of fragmenting quasiprojectiles, we wish to in- 
vestigate the possible effect of expansion of the hot 
nucleus in the temperature determination and, subse- 
quently, in the extraction of the symmetry energy coeffi- 
cient. For this purpose we employed the recently devel- 
oped expanding mononucleus model of Sobotka and Toke 
[56] . Compared with the above non-expanding Fermi gas 
framework (57.] , this model incorporates the expansion as 
a degree of freedom |6l| . The effect of the variation the 
nucleon effective mass in the level density is included in 
a manner analogous to [H^ . For a given excitation en- 
ergy E*/A, the collective expansion energy is taken into 
account. The compressibility is chosen to correspond to 
a soft equation of state. The density of the mononu- 
cleus is determined so that it maximizes the entropy. The 
temperature, in turn, is obtained from the maximum en- 
tropy state. Using this model, we calculated the values of 
the temperature for the quasiprojectiles of the reactions 
studied in this work which are summarized in Table I 
under the column T„i and later presented in Fig. 9. As 
we observe in Table I (and in Fig. 9b), the tempera- 
tures obtained with this model are systematically lower 
that those obtained from the (non-expanding) Fermi gas 
model. Despite differences in the details of the two mod- 
els, this difference may be understood qualitatively as 
the effect of expansion: in the expanding mononucleous 
model, part of the available excitation energy is allocated 
as potential energy of expansion, leaving the rest of the 
amount as thermal energy and, thus, leading to lower 
temperature. For completeness in the following discus- 
sion, both approaches for the determination of the tem- 
perature will be employed, and the respective results of 
the symmetry energy coefficient will be presented and 
discussed. 



B. Determination of the symmetry energy 
coefficient 



Having studied the isoscaling behavior and determined 
the corresponding excitation energies and temperatures, 
we currently turn our discussion to the possibility of ob- 
taining information on the symmetry energy and its de- 
pendence on excitation energy. The key quantity for 
this effort is the isoscaling parameter a. It has been 
shown [33, 112 that the isoscaling parameter a is 
directly related to the coefficient Gsym of the symme- 
try energy term E^y^n = CsymiN — Z)'^/A of the nu- 
clear binding energy. The following relation has been 
obtained in the framework of the grand canonical limit 



of the statistical multifragmentation model (SMM) y 
in the expanding-emitting source (EES) model ^fl 1^ 
and in the framework of dynamical calculations with the 
AMD model 



« = 4^((|)?-(|)^) 



(2) 



where the atomic number Z and the mass number A refer 
to the fragmenting quasiprojectiles from reactions 1 and 
2. 

In principle, Eq. [21 can serve as the basis for de- 
termining the symmetry energy coefficient Csym for ex- 
panded multifragmenting quasiprojectiles. For this pur- 
pose, the isoscaling parameter a, the isotopic compo- 
sition and the temperature of fragmenting quasiprojec- 
tiles should be determined. From the present study of 
heavy fragment isoscaling, the values of the isoscaling 
parameter a have been extracted for each reaction pair, 
as already discussed. The difference in isospin asym- 
metry between quasiprojectiles of each reaction pair is 
taken to be equal to that of the combined (fully mixed) 
system. We assumed that N/Z equilibration has been 
reached in the present reactions at the energy range 20- 
25 MeV/nucleon, as suggested by recent experimental 
studies and calculations In addition, in periph- 
eral collisions at these bombarding energies, the effect of 
pre-equilibrium emission does not appreciably change the 
difference in isospin asymmetry, as has been concluded by 
calculations using the model framework of .63j . The same 
conclusion has also been reached in the recent isoscaling 
studies of spectator fragmentation at higher energies [i^ . 
Finally, it should be noted that the effect of secondary 
decay on the values of the isoscaling parameter a in the 
excitation energy range E*/A < 3 MeV is expected to be 
small [g^. Thus, by using the experimental isoscaling pa- 
rameter a, the temperature and the difference in isospin 
asymmetry obtained as described above, we can obtain 
the values of the symmetry energy coefficient Csym for 
each case corresponding to the various values of the exci- 
tation energy, as summarized in Table I. In the following, 
we will discuss in further detail the steps involved to ex- 
tract Csym with the help of Figs. 8 and 9. 

Fig. 8 presents the experimental values of the isoscal- 
ing parameter a for the various systems studied in this 
work as a function of the difference in isospin composition 
of each pair of reactions expressed as — (^)2- The 
various symbols correspond to the reactions studied as 
explained in the caption of the figure. It is interesting to 
point out the linear relationship of the three points of the 
^^Ni data, all of which have approximately the same exci- 
tation energy of 2.9 MeV/nucleon. A similar observation 
can be made for the ^^^Xe data points, which correspond 
to a common excitation energy of 2.5 MeV/nucleon. 

Using the values of a and (^)i — (^)2j obtain the 
parameter: 



(3) 



(Z)2 
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{ZIA)\ - {Z/A)l 



FIG. 8: (Color online) Isoscaling parameter a as a function 
of {Z/A)\ — {Zj A)\ for each reaction pair studied (see text). 
Circle: *'^Kr + ^^^'^^^Sn. Square: *'^Kr + '^^■'^^Ni. Upright tri- 
angles: «*Ni + i2^'"2gj^^ 64,58j^j^ 232^j^_208p^^^ respectively in 

decreasing values of a. Inverted triangles: ^^®Xe + '^'^^'^^'^Sn, 
64,5S]\^j^ 232rpj^_i97^^^ respectively in decreasing values of a. 
The straight dotted line indicates the linear relationship of 
the ^''Ni points (see text). 



which we call "reduced" isoscaling parameter. The values 
of c with respect to the corresponding excitation energies 
are shown in Fig. 9a. Because of the common excitation 
energy for the three ^'^Ni reaction pairs (and similarly 
for the three -'^'^^Xe pairs), we plotted the corresponding 
average values of c in Fig. 9a. 

Fig. 9b shows the temperatures calculated both with 
the Fermi gas model (closed symbols) and the expand- 
ing mononucleus model (open symbols). As already 
mentioned, the temperature values calculated with the 
mononucleus expansion model are systematically lower 
(and appear to lead to a plateau with respect to excita- 
tion energy, as also discussed in [s^). 

From Eq. |21 with the definition of Eq. El for the 
symmetry energy coefficient we simply have: 
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FIG. 9: (Color online) (a) Reduced isoscaling parameter 
c = a/{{Z/A)\ - {Z/A)l) as a function E'/A. Circle: *'^Kr + 
^^-^'^^Sn. Square: «'^Kr + ^'^'^^Ki. Upright triangle: ^^Ki re- 
actions. Inverted triangle: ^''^Xe reactions, (b) Temperature 
as a function of E*/A. Open symbols: Fermi gas. Closed sym- 
bols: expanding mononucleus model ^3 ■ Symbol types as in 
(a), (c) Symmetry energy coefficient Caym as a function of 
E*/A. Open and closed symbols correspond to the two tem- 
perature estimates as in (b). Symbol types as in (a). The 
horizontal dashed line indicates the typical value of Csym for 
cold nuclei. 



Csym — ^ (4) 

The values of Csym obtained using this equation are 
shown in Fig. 9c with closed symbols (using the Fermi 
gas temperature) and open symbols (using the expanding 
mononucleus temperature). As previously discussed, the 
simple (non-expanding) Fermi gas temperatures are in 
overall agreement with the existing experimental system- 
atics, thus the extracted values of Csym should be consid- 
ered the most appropriate ones for the present analysis. 
Obviously, the values of the symmetry energy coefficient 
obtained using the expanding mononucleus temperature 
are, as expected, systematically lower. However, within 
both sets of Cgyin values, a decreasing trend with increas- 
ing excitation energy is clearly observed. In the following. 



only the values of Csym corresponding to the Fermi gas 
temperatures will be discussed further. 

From Fig. 9c (closed symbols), we observe that, at ex- 
citation energy E*/A ^ 2.0 MeV, the value of the symme- 
try energy coefficient is near (slightly lower than) the con- 
ventional value Csym,o — 25 MeV for cold (unexpanded) 
nuclei. With increasing excitation energy, however, Csym 
appears to decrease monotonically. The observed consid- 
erable decrease of the symmetry energy coefficient with 
excitation energy towards the multifragmentation regime 
is in overall qualitative agreement with the conclusions 
of Shetty et al. , Le Fevre et al. and Hcnzlova et 
al. 41]. We wish to point out, however, that the present 
study of heavy residue isoscaling reveals the gradual de- 
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crease of the symmetry energy Csym with increasing ex- 
citation energy in the range E*/A=2.0-2.9 MeV. 

The observed decrease in the fragment symmetry en- 
ergy with increasing excitation has important conse- 
quences for the formation of hot primary fragments and, 
as recently shown in [g^ , their subsequent decay. Similar 
hot nuclei are also produced in the interior of a collaps- 
ing star and subsequent supernova explosion (6Qi? ]. In 
these works it is predicted that a small decrease in the 
symmetry energy coefBcient can significantly alter the el- 
emental abundance and the synthesis of heavy elements 
in type II supernova. In light of this intimate connection, 
the present measurements and results can provide an im- 
portant testing ground to study the role of the symme- 
try energy in the formation and decay of hot fragments 
and, subsequently, to implement this knowledge in the 
simulation of the distribution of hot exotic nuclei in su- 
pernova and other hot and dense stellar environments. 
Along these lines, a thorough comparison of the isoscal- 
ing properties and the N/Z characteristics of the residues 
of the present work with the statistical multifragmenta- 
tion model (SMM) [s^l is currently underway. 



C. Final remarks and future directions 

As concluding remarks from this study, we point out 
that the present mass spectronietric data provided infor- 
mation on the isoscaling properties of heavy projectile 
fragments and, in parallel, information on the average 
excitation energy of the primary fragments (via residue 
velocity measurements). These two main experimental 
observables were used in the analysis presented in this 
work to extract the symmetry energy coefhcient at the 
corresponding excitation energy. An important advan- 
tage associated with the detailed study of heavy residues 
is that, in the respective excitation energy range E*/A 
< 3 MeV, the effect of de-excitation on the isoscaling 
parameter a is expected to be small, as discussed ear- 
lier, thus slightly affecting the procedure to obtain the 
symmetry energy coefficient. As a future experimental 
step, we wish to propose the combination of a mass sep- 
arator/spectrometer with a multidetector system capable 
of providing full acceptance and characterization for all 
fragments of the decaying projectile. Such an apparatus 
could enable the determination of the excitation energy 
on an event-by-event basis, along with temperature mea- 
surements via standard double- isotope techniques [s^ . 
This way, a detailed mapping of the excitation energy 
could be performed for various reaction pairs covering 
the whole range from low excitation energies (where only 
heavy residues are present that can be separated and 
identified with the spectrometer) to the multifragmenta- 
tion region. At each excitation energy bin, the isoscaling 
parameter a and the temperature can be obtained. In 
addition, these measurements can be supplemented by ef- 
forts to determine the de nsity (e.g., following approaches 
as discussed in [g^ and [69|'). Consequently, a correla- 



tion of the symmetry energy coefficient with excitation 
energy and, possibly, density may be obtained. Finally, 
such advanced experimental studies can be extended to 
the limits of isospin by taking advantage of present and 
future developments of rare isotope beams. 



V. SUMMARY AND CONCLUSIONS 

The isoscaling properties of isotopically resolved 
projectile residues from peripheral collisions of ^®Kr 
(25MeV/nucleon), ^"^Ni (25MeV/nucleon) and "^Xe 
(20MeV/nucleon) beams on various target pairs are em- 
ployed to probe the symmetry energy term of the nu- 
clear binding energy of hot fragments. The reactions 
of s^Kr with 64,58]xjj ^j^j i24.ii2gj^ ,^gj,g studied with 

the MARS recoil separator. The reactions of ^'^Ni and 
i^^Xe with 64,58]s^j I24,ii2gj^^ ^gjj heavier tar- 
gets (^^"^Au, 208pb, 232xh) were studied with the Su- 
perconducting Solenoid Line (BigSol) at the Cyclotron 
Institute of Texas A&M University. The present study 
focused on heavy projectile fragments produced in pe- 
ripheral and semiperipheral collisions near the onset of 
multifragment emission (E*/A = 2-3 MeV). For these 
fragments, the measured average velocities were used to 
extract excitation energies. The excitation energies, in 
turn, are employed to estimate the temperatures of the 
fragmenting quasiprojectiles within the framework of the 
Fermi gas model. The isoscaling analysis of the fragment 
yields provided the isoscaling parameters a which, com- 
bined with temperatures and isospin asymmetries, pro- 
vided the values of the symmetry energy. The extracted 
value of the symmetry energy coefficient Csym at E*/A 
~ 2 MeV is near (just below) the typical value '^25 MeV 
and is found to decrease monotonically with further in- 
crease of the excitation energy. The observed decrease of 
Csym with excitation energy is of significant importance 
to the understanding of the formation and decay of hot 
nuclei not only in nuclear multifragmentation, but also 
in supernova and other hot stellar environments. 



VI. ACKOWLEDGEMENTS 

We wish to thank A. S. Botvina for fruitful and en- 
lightening discussions at various stages of this work. We 
are thankful to L. G. Sobotka for providing us with his 
code on expanding mononucleus calculations. We are 
grateful to F. D. Becchetti and T. W. O'Donnel for valu- 
able advise and discussions during the setup and com- 
missioning of the Superconducting Solenoid Line at the 
Cyclotron Institute of Texas A&M University. Finally, 
we wish to thank the Cyclotron Institute staff for the ex- 
cellent beam quality. This work was supported in part 
by the Robert A. Welch Foundation through grant No. 
A-1266, and the Department of Energy through grant 
No. DE-FG03-93ER40773. M.V. was also supported 
through grant VEGA-2/5098/25 (Slovak Scientific Grant 



11 



Agency). 



[1] A. W. Steiner, M. Prakash, J. M. Lattimer, and P. Ellis, 

Phys. Rep. 411, 325 (2005). 
[2] C. Lee, Phys. Rep. 275, 255 (1996). 
[3] H. A. Bethe, Rev. Mod. Phys. 62, 801 (1990). 
[4] C. J. Pethick and D. G. Ravenhall, Annu. Rev. Nucl. 

Part. Sci. 45, 429 (1995). 
[5] J. M. Lattimer and M. Prakash, Phys. Rep. 333, 121 

(2000). 

[6] W. R. Hix et al., Phys. Rev. Lett. 91, 201102 (2003). 
[7] J. R. Stone et al., Phys. Rev. C 68, 034324 (2003). 
[8] A. S. Botvina and I. N. Mishustin, Phys. Lett. B 584, 
233 (2004). 

[9] K. Oyamatsu et al., Nucl. Phys. A 634, 3 (1998). 
[10] B. A. Brown, Phys. Rev. Lett. 85, 5296 (2000). 
[11] C. J. Horowitz and J. Piekarewicz, Phys. Rev. Lett. 86, 
5647 (2001). 

[12] R. J. Furnstahl, Nucl. Phys. A 706, 85 (2002). 

[13] D. T. Khoa, W. von Oertzen, and A. A. Oglobin, Nucl. 

Phys. A 602, 98 (1996). 
[14] W. Zuo, I. Bombaci, and U. Lombardo, Phys. Rev. C 60, 

024605 (1999). 

[15] M. Brack, C. Guct, and H. B. Hakansson, Phys. Rep. 

123, 276 (1985). 
[16] J. M. Pearson and R. C. Nayak, Nucl. Phys. A 668, 163 

(2000). 

[17] A. E. L. Dieperink et al., Phys. Rev. C 68, 064307 (2003). 
[18] W. D. Myers and W. J. Swiatecki, Nucl. Phys. 81, 1 

(1966). 

[19] K. Pomorski and J. Dudek, Phys. Rev. C 67, 044316 
(2003). 

[20] J. R. Stone et al., J. Phys. G, Nucl. Part. Phys. 31, 211 

(2005). 

[21] D. Lunncy, J. M. Pearson, and C. Thibault, Rev. Mod. 

Phys. 75, 1021 (2003). 
[22] D. V. Shetty, S. J. Yennello, and G. A. Souliotis, nucl- 

ex/0505011 (2005). 
[23] D. V. Shetty et al., Phys. Rev. C 70, 011601 (2004). 
[24] L. W. Chen, C. M. Ko, and B. A. Loo, Phys. Rev. Lett. 

94, 032701 (2005). 
[25] D. T. Khoa and H. S. Than, Phys. Rev. C. 71, 044601 

(2005). 

[26] B.-A. Li, C. M. Ko, and W. Bauer, Int. J. Mod. Phys. E 
7, 147 (1998). 

[27] B.-A. Li and W. U. Schroeder, eds., Isospin Physics in 
Heavy Ion Collisions at Intermediate Energies (Nova Sci- 
ence, New York, 2001). 

[28] A. Ono et al., Phys. Rev. C 68, 051601 (2003). 

[29] Y. Zhang and Z. Li, Phys. Rev. C 71, 024604 (2005). 

[30] M. B. Tsang et al., Phys. Rev. Lett. 86, 5023 (2001). 

[31] M. B. Tsang et al., Phys. Rev. C 64, 054615 (2001). 

[32] A. S. Botvina, O. V. Lozhkin, and W. Trautmann, Phys. 
Rev. C 65, 044610 (2002). 

[33] M. Veselsky, G. A. Souliotis, and M. Jandel, Phys. Rev. 
C 69, 044607 (2004). 

[34] W. A. Friedman, Phys. Rev. C 69, 031601 (2004). 



[35] G. A. Souliotis et al., Phys. Rev. C 68, 024605 (2003). 
[36] G. A. Souliotis et al., Phys. Lett. B 588, 35 (2004). 
[37] D. V. Shetty et al., Phys. Rev. C 68, 021602 (2003). 
[38] M. Veselsky, G. A. Souliotis, and S. J. YcnncUo, Phys. 

Rev. C 69, 031602 (2004). 
[39] E. Geraci et al., Nucl. Phys. A 732, 173 (2004). 
[40] A. LeFevre et al., Phys. Rev. Lett. 94, 162701 (2005). 
[41] D. Honzlova et al., nucl-cx/0507003 (2005). 
[42] A. Ono et al., Phys. Rev. C 70, 041604 (2004). 
[43] R. E. Tribble, R. H. Burch, and C. A. Gaghardi, Nucl. 

Instrum. Methods Phys. Res. A 285, 441 (1989). 
[44] W. W. Wilcke et al.. At. Data Nucl. Data Tables 25, 389 

(1980). 

[45] G. A. Souliotis et al., Phys. Rev. Lett. 91, 022701 (2003). 
[46] G. A. Souliotis ct al., Phys. Lett. B 543, 163 (2002). 
[47] G. A. Souliotis ct al., Nucl. Instrum. Methods Phys. Res. 

B 204, 166 (2003). 
[48] T. W. O'Donnell, Ph. D. Thesis, university of Michgan; 

http://www-personal.umich.edu/~twod/thesis/ (2000). 
[49] T. W. O'Donnell et al., Nucl. Instrum. Methods Phys. 

Res. A 422, 513 (1999). 
[50] G. A. Souliotis ct al.. Progress in Research, Cyclotron In- 
stitute, Texas A&M University (2003-2004), p. 11-26; ibid 

(2002-2003) p. V-5; ibib (2001-2002) p. V-19; accessible 

at: http://cyclotron.tamu.edu (2005). 
[51] H. Geissel and G. Munzenberg, Annu. Rev. Nucl. Part. 

Sci. 45, 163 (1995). 
[52] P. MoUer, J. R. Nix, and K. L. Kratz, At. Data Nucl. 

Data Tables 66, 131 (1997). 
[53] H. Madani et al., Phys. Rev. C 54, 1291 (1996). 
[54] J. Toke and W. U. Schroeder, Annu. Rev. Nucl. Part. 

Sci. 42, 401 (1992). 
[55] M. Veselsky et al., Phys. Rev. C 62, 064613 (2000). 
[56] L. G. Sobotkaet al., Phys. Rev. Lett. 93, 132702 (2004). 
[57] S. Shlomo and V. M. Kolomietz, Rep. Prog. Phys. 68, 1 

(2005). 

[58] J. B. Natowitz et al., Phys. Rev. C 65, 034618 (2002). 
[59] S. Shlomo and J. B. Natowitz, Phys. Rev. C 44, 2878 
(1991). 

[60] S. Shlomo and J. B. Natowitz, Phys. Lett. B 252, 187 
(1990). 

[61] J. Toke, J. Lu, and W. U. Schroeder, Phys. Rev. C 67, 
034609 (2003). 

[62] B.-A. Li and S. J. Yennello, Phys. Rev. C 52, 1746 (1995). 
[63] M. Veselsky Nucl. Phys. A 705, 193 (2002). 
[64] D. V. Shetty et al., Phys. Rev. C 71, 024602 (2005). 
[65] N. Buyukcizmeci, R. Ogul, and A. S. Botvina, Eur. Phys. 

Jou. A 25, 57 (2005). 
[66] A. S. Botvina and I. N. Mishustin, Phys. Rev. C 72, 

048801 (2005). 
[67] J. P. Bondorf ct al., Phys. Rep. 257, 133 (1995). 
[68] J. B. Natowitz et al., Phys. Rev. C 66, 031601 (2002). 
[69] V. Viola et al., Phys. Rev. Lett. 93, 132701 (2004). 



